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^ '. ABSTRACT 

Q I Context. The unified model of active galactic nuclei (AGN) claims that the properties of AGN depend on the viewing angle of the 

5— ( , observer with respect to a toroidal distribution of dust surrounding the nucleus. Both the mid-infrared (MIR) attenuation and contin- 

. uum luminosity are expected to be related to dust associated with the torus. Therefore, isolating the nuclear component is essential to 

Ci ' study the MIR emission of AGN. 

' ' ' Aims. This work is aimed to study the MIR emission of AGN with the highest spatial resolution available to date, isolating its con- 
tribution from extended emission. We would like to address three fundamental questions: (1) how important is the AGN contribution 

■ to the MIR spectrum? (2) does dust attenuation arise mainly from the central dusty torus associated to the AGN? and (3) how does 
^ ' spatial resolution affect these issues? 

' Methods. We have compiled all the T-ReCS spectra (Gemini observatory) available in the N-band for 22 AGN: 5 Type-1 and 17 
Type-2 AGN. The high angular resolution of the T-ReCS spectra allows us to probe physical regions of 57 pc (median). We have used 

. a novel pipeline called RedCan capable of producing flux- and wavelength-calibrated spectra for the CanariCam (GTC) and T-ReCS 
(Gemini) instruments. We have measured the fine-structure [S IV] at 10.5yum and the PAH at 1 1.3//m line strengths together with the 

' silicate absorption/emission features. We have also compiled Spitzer/IKS spectra to understand how spatial resolution influences the 

, results. We have complemented our sample with the results of 19 VISIR/VLT spectra (Paranal observatory) published by Honig et al. 

^ . (2010a) and 20 nearby highly obscured AGN (Nh > 1.5 x 10^'*cm"^) Spitzer spectra published by Goulding et al. (2012). 

' Results. The 1 1.3/jm PAH feature is only clearly detected in the nuclear spectra of two AGN, while it is more common in the Spitzer 

• • ' data. For those two objects the AGN emission in NGC 7130 accounts for more than 80% of the MIR continuum at 12;um while in the 

^ , case of NGC 1808 the AGN is not dominating the MIR emission. This is confirmed by the correlation between the MIR and X-ray 

■ continuum luminosities. The [S IV] emission line at 10.5yum, which is believed to originate in the narrow line region, is detected in 



most AGN. We have found an enhancement of the optical depth at 9.1 jim (T9 7) in the high-angular resolution data for higher values 
of Nh. Clumpy torus models reproduce the observed values only if the host-galaxy properties are taken into account. 
Conclusions. 

Key words. Accretion, accretion disks - Galaxies: active - Galaxies:Seyfert - Infrared:galaxies - Techniques:imaging spectroscopy 



1. Introduction 

The radiation from active galactic nuclei (AGN) is due to accre- 
tion onto a super-massive black hole (SMBH). The optical spec- 
tra of AGN have been classified into two main classes according 
to the existence (Type- 1 ) or not (Type-2) of broad permitted lines 
(FWHM> 2000 km 5"'). The so-called unified model proposes 
that both types of AGN are essentially the same objects viewed at 
different angles (Antonucci, 1993; Urry & Padovani, 1995). An 
optically-thick dusty torus surrounding the central source would 
be responsible for blocking the region where these broad emis- 
sion lines are produced (the so-called broad line region, BLR). 



* Juan de la Cierva Fellow (e-mail: omairagm@iac . es) 



Therefore, Type-2 AGN are essentially Type-1 AGN blocked by 
the dusty torus along the line of sight (LOS) to the observer. 
More elaborated models suggest that the same observational 
constraints can also be explained with a torus made of discrete 
dense molecular clouds following a certain radial distribution 
(EHtzur, 2008; Honig et al., 2006; Honig & Kishimoto, 2010b; 
Nenkova et al., 2008a,b; Schartmann et al., 2008). 

Mid-infrared (MIR) spectroscopy is a powerful tool to ex- 
amine the nature of AGN, as well as the contribution of star- 
formation activity. Weedman et al. (2005) found a variety of 
spectral features displayed in the MIR spectra of eight classi- 
cal AGN with Spitzer data (see also Rodriguez-Espinosa et al. , 
1987, 1997). Large samples of AGN observed with Spitzer have 



1 



O.Gonzalez-Martin et al.: Dust in AGN 



Num 


Name 


D(Mpc) 


Slit width(arcsec)/Size(pc) 


Type 


Optical features 


log(NH) 


log(L(2-10keV)) 


X-ray Ref 


b/a 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


1 


NGC 1365 


17.9 


0.35/30 


S1.5 




23.6 


42.5 


(1,11) 


0.5 


2 


NGC 1386 


16.2 


0.31/24 


S2 


E 


>25.0 


42.6 


(2,3) 


04 


3 


NGC 1808 


11.5 


0.35/19 


S2 




22.9 


40.4 


(3,12) 


0.6 


4 


NGC 3081 


32.5 


0.65/102 


S2 




23.9 


42.5 


(2) 


0.8 


5 


NGC 3094 


44.7 


0.35/76 


S2 










0.7 


6 


NGC 3281 


45.7 


0.35/77 


S2 


D,E 


24.3 


43.4 


(1,12) 


0.4 


7 


NGC 44 18 


31.1 


0.35/53 


S2 


E 


>25.0 


41.2 


(4) 


0.4 


8 


NGC 4945 


4.0 


0.65/13 


S2 


D,E,M 


24.7 


42.3 


(5,12) 


0.4 


9 


Centaurus A 


3.7 


0.65/12 


S2 


D,E 


23.7 


41.9 


(1,13) 


0.4 


10 


NGC 5135 


58.6 


0.70/199 


S2 




>25.0 


43.1 


(2) 


0.7 


11 


Circinus 


4.2 


0.35/7 


S2 


M 


24.6 


42.6 


(6,2) 


0.5 


12 


NGC 5506 


28.8 


0.35/49 


S1.9 


D,E 


22.5 


43.0 


(1,3) 


0.2 


13 


NGC 5643 


16.9 


0.35/29 


S2 




23.8 


42.6 


(7,12) 


0.9 


14 


NGC 5728 


29.1 


0.35/49 


S1.9 




23.6 


43.3 


(6) 


0.6 


15 


IC4518W 


69.6 


0.70/236 


S2 


E,M 


24.3 


42.6 


(8,13) 


0.1 


16 


hi>0 iU3-G35 


56.9 


0.35/96 


S2 


E 


23.2 


43.4 


(1,12) 


0.4 


17 


IC5063 


48.6 


0.65/153 


S2 


D 


23.3 


43.3 


(1,12) 


0.7 


18 


NGC7130 


69.2 


0.70/235 


L2 


M 


>25.0 


43.1 


(10,2) 


0.9 


19 


NGC7172 


33.9 


0.35/57 


S2 


D 


22.9 


42.8 


(1,3) 


0.6 


20 


3C445 


239.3 


0.35/406 


S1.5 




23.5 


43.9 


(3) 


0.9 


21 


NGC 7479 


33.8 


0.35/57 


L1.9 


D,M 


24.3 


42.0 


(3) 


0.8 


22 


NGC 7582 


20.6 


0.70/70 


S2 


D,E 


23.0 


42.6 


(2) 


0.4 



Table 1. Sample properties. The size of the region observed by T-ReCS is computed using the slit width. Optical features: (D) Dust features across the nuclear region; (E) Edge-on galaxies; 
(M) Merger. The references for the X-ray results are: (1) Tueller et al. (2008), (2) Marinucci et al. (2012), (3) Brightman & Nandra (201 1), (4) Maiolino et al. (2003), (5) Itoh et al. (2008), 
(6) Goulding et al. (2012), (7) Guainazzi et al. (2004), (8) Pereoa-Santaella et al. (2011), (9) Moran et al. (1996), (10) MuUaney et al. (2011), (11) Risaliti et al. (2009), (12) Comastri 
(2004), (13) de La Rosa et al. (201 1), and 



been studied in detail to quantify the properties of AGN at MIR 
(e.g. Deo et al., 2007; Hao et al., 2007; Sturm et al., 2006). 

Among the most noteworthy results found at MIR wave- 
lengths, Wu et al. (2009) found that the silicate features at 9.1 
and 18/im in Type-1 Seyfert galaxies are rather weak, while 
Type-2 Seyfert galaxies are likely to correspond to strong sil- 
icate absorptions. Furthermore, the deepest sihcate absorption 
features show the highest values for the hydrogen column den- 
sity (Nh) in X-rays (Shi et al., 2006). The silicate emission 
feature in Type-1 AGN is much weaker than expected from 
smooth dusty tori. The dumpiness of the dusty torus explains 
better the overall suppression of this feature (Nenkova et al., 
2002, 2008a). However, clumpy models do not predict strong 
silicate absorption features because illuminated and shadowed 
clouds contribute to the observed spectrum. In fact, the sili- 
cate feature depth typically seen in some Compton-thick AGN 
(Nh > 1.5 X 10^'^cm"^) observed with Spitzer is similar to that 
observed in Type-1 AGN (Goulding et al., 2012). However, the 
nature of these silicate feature is not clear since the host-galaxy 
could be another contributor to the dust attenuation seen in the 
MIR (see also Alonso-Herrero et al., 2011; Deo et al., 2007). 

Another prominent feature noticeable in the MIR is the poly- 
cyclic aromatic hydrocarbon (PAH) emission. This is usually 
prominent in starburst galaxies, with a good correlation between 
the strength of the PAHs and the IR luminosity in starbursts 
(Brandl et al., 2006). This correlation appears however to be 
weak or absent in AGN (Weedman et al., 2005). PAHs form 
in the photodissociation regions (PDRs) in star forming clouds, 
but these PDRs can easily be destroyed by the intense radiation 
field from the AGN, even if it is a vigorous starburst (Genzel 
et al., 1998; Laurent et al., 2000). Goulding et al. (2012), using 
low resolution Spitzer spectra, found that Compton-thick AGN 
show strong PAH features, suggesting that the PDRs can survive 
the AGN radiation field if properly shielded. Diamond-Stanic & 
Rieke (2012) found a relation between the accretion rate of AGN 
and the SFR derived from PAHs. They suggest that the gas fuel- 
ing the SMBH is also used for the star-formation. 

However, the relatively low-spatial resolution of Spitzerf[RS> 
spectra (~ 3 arcsec) include a non-negligible contribution of the 



host galaxy and/or of circumnuclear star formation. Thus, a good 
isolation of the torus emission is needed to understand the pure 
AGN contribution to the MIR spectrum. The best isolation possi- 
ble to date can be obtained using IR interferometry (e.g. Tristram 
et al., 2007, 2009). However, this can only be done for a few 
bright AGN. As for the rest of AGN the medium- to large-size 
ground-based telescopes oflfer an opportunity to better isolate the 
contribution of the AGN within a few tens of parsecs. Some 
individual works have been presented showing the importance 
of spatial resolution to isolate the AGN from the host galaxy 
(e.g. NGC 1068, NGC 1365, NGC 3281 and the Circinus galaxy; 
Alonso-Herrero et al., 2012; Mason et al., 2006; Roche et al., 
2006; Sales et al., 2011, respectively). A large sample was pre- 
sented by Honig et al. (2010a), including 19 AGN observed with 
the VLT/VISIR MIR spectrograph with spatial resolution better 
than ~100 pc. They found moderately deep silicate absorption 
features. However, this sample is free of objects with high in- 
clination angles or dust lanes, and include only three Compton- 
thick AGN. The MIR luminosity seems to be quite isotropic, 
correlating with other AGN luminosity tracers (see also Gandhi 
et al., 2009; Levenson et al., 2009; Mason et al., 2012; Ramos 
Almeida et al., 2007). 

In this paper, we compile all the AGN observed by T-ReCS at 
the Gemini observatory. The sample comprises 22 AGN, includ- 
ing 9 Compton-thick sources. Our main aims are: (1) to charac- 
terize the AGN pure emission at MIR frequencies; (2) to under- 
stand the nature of the attenuation seen at MIR frequencies as 
traced by the optical depth at 9.7 fim; (3) to compare T-ReCS and 
Spitzer spectra to understand the limitations introduced by the 
lower spatial resolution in Spitzer spectra. 

For the data processing and analysis we have used a newly 
developed pipeline called RedCan^ This pipeline is able to 
produce flux-calibrated images and ID spectra for both the T- 
ReCS (De Buizer & Fisher, 2005; Telesco et al, 1998) and the 
CanariCam^ (Telesco et al., 2003) instruments without user in- 
teraction. This pipeline has been produced to facilitate a better 

' The CanariCam data reduction pipeline. 

^ CanariCam is a mid-infrared (7.5 - 25/jm) imager with spectro- 
scopic, coronagraphic, and polarimetric capabilities, which is mounted 
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Fig. 1. Flowchart of the RedCan procedure. Diamonds are the steps of the imaging procedure while polygons indicate the steps of the spectroscopic 
procedure. The empty flag is the input level and the dark gray flags indicate the outputs levels. The dashed arrow is set to remark that RedCan 
continues from imaging to spectroscopic procedures only when spectra are included in the initial dataset. 



exploitation of the CanariCam data and can be freely accessed 
(see Section 3). 

2. Sample selection 

We have searched for all the T-ReCS (Gemini South observa- 
tory) AGN spectra available in the Gemini public archive^ until 
June 2012. The sample comprises 22 AGN observed as part of 
~10 observational projects. All of them are observed in the N- 
band. Note that, since T-ReCS is not expected to be available 
after 201 3 A this sample includes most of the AGN ever to be 
observed with this instrument'*. 

Table 1 shows the main observational details of the sample. 
Seventeen objects are Type-2 and five are Type-1 AGN. The slit 
width used for the spectroscopy results in a spatial resolution be- 
tween ~ 7-236 pc for all the objects except for 3C 445 (406 pc), 
which is by far the most distant object in the sample. This is 
the largest sample of high-spatial resolution spectra ever put to- 
gether, especially when including the 19 AGN (VISIRA'LT) re- 
ported by Honig et al. (2010a). Our sample includes 9 Compton- 
thick AGN (see Nh in Table 1), making it complementary to the 
AGN sample in Honig et al. (2010a) (where only 3 Compton- 
thick AGN are included). Note, however, that our sample is un- 
biased though not complete in any sense since the observations 
come from several different proposals with differing goals. It is 
worth noticing that this sample appears to contain a high fraction 
of sources with large silicate absorption features (see Section 6). 

3. Data reduction: the CanariCam data reduction 
pipeline (RedCan) 

Data reduction pipelines are an important aspect of modem as- 
tronomy. Indeed large telescope instruments are becoming fairly 
complex, hence the prompt availability of science data is cru- 
cial for securing fast scientific turnout from today's large and 
expensive facilities. Unfortunately, this is not always the case, 

at one of the Nasmyth foci of the GTC telescope at El Roque observa- 
tory. 

' http://www3.cadc-ccda.hia-iha.nrc- 
cnrc.gc.ca/cadcbin/gsa/wdbi.cgi/gsa/gsa _science/form 
Except for three objects that are not public yet. 



even in successful facilities. The aim of our group is to produce 
a new and advanced pipeline for MIR imaging and spectroscopic 
data. In the forthcoming years large amounts of GTC/CanariCam 
data will be acquired by our group. This includes data from both 
guaranteed CanariCam time (~100 hours, PI C. Telesco) and 
the ESO/GTC proposal led by A. Alonso-Herrero (~180 hours, 
ESO/GTC 182.B.2005). It is therefore important to have a data 
reduction pipeline that facilitates the use of the data as soon as 
possible after being acquired at the telescope. This is the purpose 
of RedCan, the Canaricam data reduction pipeline we have de- 
veloped. RedCan is a well tested reduction and calibration pro- 
cedure to ensure reliable data products from CanariCam. 

Thanks to the similarities between the T-ReCS and 
CanariCam instruments, RedCan can be used for data from both 
instruments in both imaging and spectroscopy modes (Q and N 
bands in medium resolution). 

RedCan has been written combining the Gemini IRAF pack- 
age^) within PyRAF, together with ad hoc IDL routines devel- 
oped for this purpose. CSHELL language is used for the main 
routine of the pipeline to connect IDL and Python routines. ^. 
Fig. 1 shows a chart flow of the reduction steps done by RedCan 
to produce the final reduced data. RedCan is based on the pro- 
cedure described by Dfaz-Santos et al. (2010) (see also Dfaz- 
Santos et al., 2008) and the main input file is an ASCII file with 
a list of observations. A full description of the pipeline procedure 
is given below. 

3. 1 . Identification of files 

A full data set from a target observation consists of: (1) telluric 
standard observations, (2) target observation, and (3) acquisi- 
tions of both standard and target images. Standards and targets 



^ The released version of the Gemini IRAF package is an 
external package layered upon IRAF and is available to users 
and other interested parties ( http://www.gemini.edu/sciops/data-and- 
results/processing-software 

' The RedCan package can be retrieved here: 

https://www.dropbox.com/sti/w4z8buo2ewrhhvj/Kklf2RYVe8/RedCan.tar.gz. Note 

that the pipeline is of free access although it will not be supported. A 
manual is also available that includes documentation on the software 
requirements and installation procedures. 
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can be in imaging and/or spectroscopic modes. Moreover, the ac- 
quisition images in the case of spectroscopic data can be either 
imaging or imaging throughout the slit. 

RedCan uses keywords in the fits headers to identify the type 
of each observing block, which are used accordingly throughout 
the pipeline. In particular the keyword obsclass together with the 
instrumental configuration keywords (i.e. slit, grating, filterI, 
filter2, etc.) and observational keywords (i.e. ontime, ra, dec, 
UT, date) are used to identify the following information: 

- Acquisition image associated to a target/standard. RedCan 
sorts the files according to the observed time and selects 
as acquisition images those immediately before the tar- 
get/standard. Acquisition images of the objects through the 
slit are used to compute slit losses for the spectroscopic mode 
(see Section 3.8). These images are also flux calibrated dur- 
ing the process so they can be used for science purposes. 

- Standard star image/spectrum associated to a target. RedCan 
selects within the same date of observation the standard im- 
age/spectrum observed with the same configuration used for 
the target (i.e. same filter and grating configuration with the 
same slit size). Note that if more than one observation of 
standard stars fulfill this criterion, the one closer in time is 
chosen. This association is used to flux calibrate both images 
and spectra (see Sections 3.4 and 3.8). 

3.2. Flat fielding 

The RedCan pipeline allows to subtract a master flat computed 

as an average of the observations of the sky for each chop in- 
tegration. This master-flat is computed using the tasks tprepare 
and MisTACK within the midir/IRAF package while the flat sub- 
traction is implemented in IDL. Diaz-Santos et al. (2008) con- 
cluded that flat-fielding can improve the quality of the final im- 
ages. However, the flat-fielding procedure introduces noise into 
the final image (see e.g. the Gemini web page^). We have de- 
tected that using flat-fielding in the spectroscopic mode can lead 
to foUow a wrong trace of the spectra due to the inclusion of pat- 
terns at the borders of the array. This step is included in RedCan 
as an option for the user, although the default is set to skip it. We 
have skipped this option in the data presented in this paper. 

3.3. Stacking 

Accurate cancellation of the sky and telescope background is 
most readily achieved through the technique of chopping and 
nodding. The sky is very bright and highly variable at MIR fre- 
quencies when observed from the earth. To get rid of the sky, 
the chop technique observes frequently a nearby sky position 
(ofl^-source chop) moving the secondary mirror of the telescope 
to subtract from the target position (on-source chop). However, 
this chopping technique leaves a background resulting from the 
different paths followed by the light depending on whether the 
telescope is pointing to the on-source or the ofl-source positions. 
The nodding technique places the source at the nod positions 
moving the actual telescope axes to remove this residual. Thus, 
the chop-nod technique allows the subtraction of the sky back- 
ground as well as the background from the telescope itself. This 
technique also removes the array 1/f noise. The data are then 
stored automatically in save-sets. Prior to stacking, RedCan au- 
tomatically rejects chop-nod pairs of images affected by noise 
patterns or glitches due to instrumental or background issues 

^ http://www.gemini.edu/?q=node/ 101 39#flat 



Name 


ObsID 


Expt. 


DatE 


P.A, 


STD* 


FWHM 
(8.7/im) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


NGC 1365 


0061 


624 


090829 


15 


26967 


0.14 


NGC 1386 


0210,12 


2x606 


060929 





23319 


0.29 


NGC 1808 


0076 


981 


090829 


45 


26967 


0.20 


NGC 3081 


0094 


606 


060125 





81797 


0.20 




0113 


606 


060304 


350 


93813 


0.29 


NGC 3094 


0075,76 


2x309 


040506 





81797 


1.10 


NGC 3281 


0120 


981 


090328 


315 


HR3438* 


0.17 




0033,35 


2x981 


090406 


315 


146791 


0.20 


NGC 4418 


0041 


309 


040507 


40 


128620 


0.25 




0046,47 


2x309 


040508 


40 


128620 


0.13 




0025-28 


4x327 


080628 


30 


99998 


0.30 


NGC 4945 


0055,56 


2x303 


060417 


45 


110458 


0.26 


Centaurus A 


0198 


517 


040131 





123139 


0.22 


NGC 5 135 


0098 


633 


060306 


30 


101666 


0.19 




0106,07 


2x633 


060310 


30 


101666 


0.16 


CIRCINUS 


0084,85 


2x309 


040506 


100 


81797 


1.10 




0095-97 


3x309 


040506 


10 


81797 


1.10 


NGC 5506 


0107,08 


2x309 


040506 





81797 


1.10 


NGC 5643 


0057-60 


4x303 


070605 


80 


123139 


0.21 




0057-59 


3x327 


080615 


80 


123139 


0.50 




0060-62 


6x327 


080615 


80 


139127 


0.40 


NGC 5728 


0046,47 


2x921 


090406 


215 


146791 


0.18 


IC4518W 


0081-83 


3x633 


060417 


5 


123123 


0.22 


ESQ 103-G35 


0065-67 


3x318 


080616 


30 


156277 


0.24 




0068-70 


3x318 


080616 


30 


171759 


0.51 


IC5063 


0078 


606 


05-05-21 


303 


81797 


0.20 




0046 


606 


050707 


303 


81797 


0.20 




0047 


303 


050707 


303 


81797 


0.20 




0048 


302 


050707 


303 


81797 


0.20 




0092 


302 


050709 


303 


81797 


0.20 


NGC 7130 


0061 


633 


050918 


348 


190056 


0.30 




0062 


633 


050918 


348 


219784 


0.35 




0111-113 


3x633 


060704 


348 


190056 


0.14 




0057 


633 


060829 


348 


190056 


0.24 




0107,08 


2x633 


060916 


348 


190056 


0.19 




0097,98 


2x633 


060925 


348 


190056 


0.25 




0099 


633 


060925 


348 


219784 


0.29 


NGC 7172 


0120-22 


3x309 


040506 


60 


81797 


1.10 




0077-80 


4x327 


080615 


90 


200914 


0.29 




0041,42 


2x318 


080625 


90 


200914 


0.34 




0063-66 


4x327 


080627 


90 


200914 


0.25 


3C445 


0098,99 


2x921 


100606 


300 


206445 


0.24 


NGC 7479 


0055-57 


3x327 


080629 


10 


206445 


0.41 




0058-60 


3x327 


080629 


10 


1255 


0.23 


NGC 7582 


0174,75 


2x606 


050718 





BS8556* 


0.20 



Table 2. Observational details of the T-ReCS spectra. *STD: Standard star used to per- 
form flux-calibration. All the standard stars are named with their HD sequence num- 
ber except those marked with a star. The ObsID (Col. 2) includes the last 4 digits of 
the name of the observing file. The first 8 digits are the date of the observation (i.e. 
"S"H-ObsID4-"S20"+YYMMDD). 



during the observations. This rejection is done by default exclud- 
ing any chop-nod pair above or below 5-cr of the mean value 
of the subtracted chop-nod pair of images (task tbackground 
within Gemini/midir). However, this can be executed manually 
and the threshold level can be changed if desired. Both the on- 
source and the ofi'-source chop images are averaged together 
(task miprepare) and each nod pair combined (task mistack). The 
resulting image shows certain types of pattern noise as horizon- 
tal or vertical stripes. Those that are constant along the image are 
removed with the task miclean within the gemini/midir package. 
Finally each image is divided by the frame time so the final im- 
age is in units of "adu/s". 

The final product is a single data cube per observation with 
three planes: on source, off source, and on-off source (i.e. back- 
ground subtracted) images. This procedure is done for every im- 
age (i.e. acquisitions, standards, and targets). Note that the sec- 
ond plane is needed for wavelength caUbration in the spectro- 
scopic mode (see Section 3.5). 
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3.4. Imaging flux calibration 

Images resulting from the stacking process described above will 
now be flux calibrated. This flux calibration is performed for 
both acquisition and target images using associated telluric stan- 
dards stars. 

RedCan contains the catalog of synthetic flux-caUbrated 
spectra of standards stars published by Cohen et al. (1999). An 
identification of the standard star is done using the coordinates 
given in the header of the observation. When the standard star is 
not part of the catalog of Cohen, a reference standard star is used 
and a waming is shown in the screen and recorded in the out- 
put tables and log of the reduction process. Then, the synthetic 
spectrum of the star is convolved with the telescope through- 
put and filter bandpass. Aperture photometry of the standard 
star is performed and a new keyword called fcalstd is placed 
in the standard star header which provides the ratio between 
the number of counts per second and the real flux of the star. 
Then aperture photometry techniques are applied to the target 
images and a table with the final fluxes for the set of images is 
recorded. Note that aperture radii of the standard star and target 
are fixed to a radius of 2.7 and 0.9 arcsec respectively. These 
radii include more than 95% of the flux at the diflraction limit in 
both T-ReCS (FWHM~0.35-0.37 arcsec ^) and CanariCam data 
(FWHM~0. 26-0.42 arcsec and above 90% with good seeing 
conditions (optical seeing below 1 arcsec). The background is 
extracted using annuli with intemal and extemal radii of [4.5, 
5.3] and [1.8, 2.25] arcsec, respectively. Finally, flux calibrated 
images are also produced for target and standard stars using the 
conversion factor stored in fcalstd. 

This method is flux accurate to 10% (Dfaz-Santos et al., 
2008). However, we note that this aperture photometry is ac- 
curate enough for point-like sources while a more sophisticated 
methodology is needed for extended sources (see Section 4.1). 
RedCan finishes here if the input Ust of files does not include 
any spectroscopic data. It continues to the next steps otherwise. 

3.5. Wavelength calibration 

RedCan uses sky emission features in the reference stacked spec- 
tra to wavelength calibrate all the spectra (standard stars and tar- 
gets). First, the sky spectrum is stacked along the spatial axis. 
The pipeline looks for the local maxima in this stacked spectrum 
as local peaks within the 10 closest pixels. A gaussian profile is 
fitted to the 1 1 pixels centred in the local maxima to accurately 
determine the position of the centre. This method is rather sen- 
sitive to broad features such as those from the MIR sky. 

The next step is to identify these local maxima with 
known features. For that purpose the pipeline contains a set of 
wavelength-calibrated sky templates for the N and Q band and 
for T-ReCS and CanariCam. These observational templates are 
obtained from the best weather conditions data from the Gemini 
archive and commissioning data, for T-ReCS and CanariCam 
data, respectively. After selecting the associated template ac- 
cording to the grating and instrument, a grid of offsets and dis- 
tortions is applied to it. The best offset and distortion for the 
sky-spectrum is chosen by computing the minimum x^- Using 
the wavelength calibration of the templates, RedCan identi- 
fies as many features as possible in the sky-spectrum. A re- 
identification of these lines along the spatial axis is made each 



* http://www.gemini.edu/?q=node/10582#trecspsf 
' http://www.gtc. iac.es/en/pages/instrumentation/canaricam/data- 
commissioning.php #Image_Quality 



10 pixels. A linear fit is then applied along the wavelength direc- 
tion and third and fourth order fit are appUed along the spatial 
direction for T-ReCS and CanariCam data respectively. The fi- 
nal calibration is stored in a 2D wavelength calibration array. 

The N-band spectral resolution is R=100 and R=175 for T- 
ReCS and CanariCam respectively. The accuracy of the wave- 
length calibration is always much better than the spectral resolu- 
tion of the instruments. 

3.6. Trace determination 

Before extracting the spectra the traces of standard stars and tar- 
gets are computed. The pipeline stacks the spectrum along the 
wavelength axis every 10 pixels and determines the centre and 
width of the trace by fitting it to a MOFFAT profile (the mpfit- 
peak'" function of IDL). The resulting widths of the trace along 
the spatial direction are fitted to a polynomial. A first and second 
order degree polynomials are tested against the N-band spectra 
to include the possibiUty that the data might not be diffraction- 
limited at bluer wavelengths. For the width and center of the 
trace in the Q-band a linear fit is used for both T-ReCS and 
CanariCam. The centre and width of the trace for each obser- 
vation are saved in its header with the keywords trctra, trctrb, 

TRFWHMA, TRFWHMB, and TRFWHMC. 

3. 7. Spectral extraction 

RedCan extracts the spectrum by using the 2-D final spectrum of 
the target and standard star and the wavelength calibration ma- 
trix (see Section 3.5). For the centre of the extraction region the 
pipeline uses the centre of the trace for each object recorded in 
the header (see Section 3.6). For the width of the extraction re- 
gion two considerations can be made: (1) extended or (2) point- 
hke sources^'. 

For point-like sources the width along the spatial direction 
as a function of wavelength is set to vary as the trace of the as- 
sociated standard star. The user can specify the number of times 
the width of the standard to be used. To extract the unresolved 
component this number should be set to 1 . For extended sources 
the width does not vary with wavelength and the user can spec- 
ify the number of pixels to be extracted. Additionally, a number 
of offsets (in pixels) from the centre of the target (i.e. centre of 
the trace) can be specified to extract one or several regions out 
of the centre of the target. The pipeline can extract either several 
widths for an offset or several offsets for a given width. 

For our AGN sample we have selected the point-Uke source 
option for unresolved sources using the width of the standard star 
to isolate the nuclear component from the extended emission. 

3.8. Spectral flux calibration and combination of spectra 

Each spectrum is flux-calibrated using the ratio between the ob- 
served associated standard and its corresponding flux calibrated 
template given by Cohen et al. (1999). For the extraction of the 
standard a radius of 4 times the sigma of the trace of the standard 
is used to ensure the collection of the whole flux of the standard. 

For extended sources the final flux-calibrated spectrum is ob- 
tained by dividing the spectrum (in counts) by this ratio. For 
point-Uke sources, in order to obtain a well flux-calibrated spec- 
trum, losses due to the width of the slit (slit-losses) and due to 
the aperture used (aperture-losses) are taken into account. The 

http://cow.physics.wisc.edu/ craigm/idl/idl.html 
" These two options can be selected by the user. 
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slit losses are automatically computed using the acquisition im- 
ages for the associated standard star. Two images must be present 
within the observations, before and after including the slit in the 
optical system. The pipeline computes aperture photometry in 
these two images and put the ratio of the counts of the obser- 
vations without slit divided by the counts using the slit in the 
header adding the new keyword slitcorr in the observation cor- 
responding to the spectrum of the star. When the user selects the 
option point-like the final spectrum is multiplied by this value 
to obtain the final flux-calibrated spectrum. Note also that slit- 
losses are wavelength-dependent for point-like sources because 
the angular size of the PSF in a diffraction limited observation 
increases linearly with wavelength. Therefore the loss ratio cal- 
culated using the acquisition images, with and without slit, is 
only a first order correction, since it is based on the averaged flux 
of the source as measured through the filter used for each acqui- 
sition and thus is a constant value. While this is not important 
for point-source extractions, it is a source of uncertainty for ex- 
tended sources. When observing unresolved sources, both stan- 
dard and target are aff'ected by the same wavelength-dependent 
slit loss, which is automatically accounted for when the spec- 
trum of the target is divided by that of the standard. On the other 
hand, extended sources are assumed as not having any loss due 
to either the slit or aperture. Therefore the standard star intro- 
duces an uncertainty in the shape of the extracted spectrum of 
the target even when the standard is corrected by the "constant" 
loss ratio. 

For the aperture losses the standard star is extracted with the 
same width as the target. These aperture losses depend on the 
wavelength and they are computed as the ratio between the stan- 
dard spectrum (with 4 times the standard trace) and the spectrum 
extracted with the same trace as the target. 

The errors are estimated as the sum of the statistical error 
(squared root of the number of counts) plus a 15% systematic 
errors due to the flux calibration procedure (see Diaz-Santos et 
al., 2008). 

The observations are grouped together into targets using the 
keyword objname in the headers. All the flux- and wavelength- 
calibrated spectra of the same target are combined into one. This 
final spectrum is given at rest-frame by using the redshift ob- 
tained from NED'^. For this the name of the target is used. When 
the object is not found in NED (either because it is not extra- 
galactic or the name is not well spelled in the header), a redshift 
equal to zero is used. Note that the applied redshift is stored in 
an ascii file. 

Table 2 shows the observational details for our sample. In 
particular the information on the standard star used to flux- 
calibrate each observations and the number of observations com- 
bined to produce the final spectrum of each target. 

4. Data analysis 

4. 1 . MIR images 

RedCan produces flux-calibrated images (see Section 3.4) using 
the associated standard star images and their theoretical spectra 
given by Cohen et al. (1999). Fig. 2 (left) shows the resulting 
flux-calibrated image for NGC 3094. Most of the N-band im- 
ages are acquisition images associated to the target. For that 
reason most of the sources have imaging in the N-band (ex- 
cept Centaurus A, NGC 5506, and NGC 5728). Q-band filter im- 
ages are available for 5 objects while the narrow band filters Si2 



http ://ned. ipac . caltech. edu/ 
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Fig. 2. Si6 image at 12.3fim of NGC 3094 before (left panel) and after 
(right panel) the PSF subtraction. In the latter, one can clearly see the 
nuclear emission in red. The dashed contour is the best MOFFAT fit for 
the standard star. The color scale is in units of ly. 



(8.8fim) and Si6 (I2.3fim) are available for 9 and 3 sources, re- 
spectively. All sources have an image at least in one of the filters. 

Our aim is to study the central region morphology of the 
sample in the MIR and isolate the AGN contribution. We have 
used the closest (in time) standard star observed with the same 
filter as representative of the PSF of a point-like source. This star 
is fitted to a 2D MOFFAT profile (using the "mpfit2dpeak" IDL 
function). We have scaled this PSF function to the peak of the 
nuclear source and subtracted the nuclear component. Note that 
this method would overestimate the nuclear flux of the source if 
extended. However, it is not our aim to measure fluxes (they are 
estimated with the spectra) and this approximation is robust at 
determining the number of extended sources. 

The right panel in Fig. 2 shows the residual after the sub- 
traction of the PSF (dashed contours) for NGC 3094. We con- 
sider a source as extended if the image after the PSF subtraction 
shows an integrated flux higher than 3 times the standard devia- 
tion of the background. Column 2 in Table 3 indicates whether 
the sources are extended (E) or point-like (P). 

In general, for the sources in our sample, the central mor- 
phology is extended. Only 6 AGN show point-like nuclear emis- 
sion, namely NGC 1365, Centaurus A, NGC5728, IC4518W, 
ESO 103-G35, and NGC 7479. However, Alonso-Herrero et al. 
(2012) found an extended morphology for NGC 1365 using N- 
band images with T-ReCS. Their proprietary data have a higher 
signal-to-noise since our image is only the acquisition image for 
the spectrum. Thus, the S/N influences the identification of the 
extended emission in NGC 1365. In the cases where the source 
is not extended the nuclear emission comes from less than 7-236 
pc, depending on the distance to the galaxies. The sources ob- 
served using several filters show similar morphologies in all the 
filters. NGC 4945 is the only AGN without an unresolved point- 
like source. Overall, we find that the morphology of the nuclear 
source has not clear relation with the AGN type. 

4.2. MIR spectra 

Fig. 3 shows the nuclear spectrum of NGC 7582 (open dots and 
gray area). The spectra of the entire sample are included in ap- 
pendix A. The shaded regions show the systematic errors of the 
extracted spectra. The spectra show systematic eiTors of less than 
10% of the flux in all the objects except for NGC 5728, in which 
the error is much higher (~ 40%). We have excluded NGC 5728 
from this spectral analysis due to the lack of good quality data. 
Moreover, NGC 4945 is also excluded due to the lack of a nu- 
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Table 3. Main results. T9.7 is measured in the 8-13/jm band except when using PAHFIT in the Spitzer spectra for which we have used the full 5-35/jm band. (P) Point-like source and (E) 
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Fig. 3. T-ReCS (open dots) and Spitzer/IRS (red-solid line) spectra of 
NGC 7582. The errors are shown as a grey area. Interpolated continuum 
and estimated silicate feature depth are shown as a black-solid (red- 
dotted) line and vertical solid blue (red dotted-dashed) line for the T- 
ReCS (Spitzer/IRS) spectrum. The appendix A includes the spectra of 
the sample. 



clear component. Thus, our sample of nuclear spectra comprises 
20 AGN. Our sample has three objects in common with the 
sample reported by Honig et al. (2010a) using VISIR/VLT data 
(NGC 5643, IC 5063, and NGC 7582). The shape and detected 
features in the three VISIR spectra in common with our sample 
are consistent with those reported with T-ReCS. The total flux 
calibration is consistent with ours except for IC 5063, where we 
find a flux level twice fainter than that reported by Honig. 



The low spatial resolution spectra (Short-Low module) 
Spitzer spectra in the same spectral range as T-ReCS are also 
shown in the appendix when available (red thick line). Spitzer 
data were retrieved from the Cornell Atlas of SpitzerfTRS 
Sources (CASSIS Lebouteiller et al., 2011). Note that these 
spectra are extracted using the optimal extraction region to 
ensure the best signal-to-noise ratio. IC4518W was reduced 
with the Smart package'"' (version 8.2.6 Higdon et al., 2004; 
Lebouteiller et al., 2010) since it was not available through 
CASSIS. There were Spitzer/IRS spectra for all the AGN in our 
sample except for NGC 3094. 

The observed apparent optical depth at A ~ 9Jfj.m can 
be measured as the ratio between the expected continuum 
flux (fgj^com) and the observed flux in the feature (fgj^obs, see 
Georgantopoulos et al., 2011; Levenson et al., 2007; Shi et al., 
2006; Spoon et al., 2007): 

T9.7 = ln(f()j^cont/f9J,obs) (1) 

Sources with absorption (emission) features show positive 
(negative) optical depths. To estimate f^.j^cont we have fitted the 
continuum to a linear relation to the %.5pm and 12.5jum flux for 
each object. Fig. 3 shows the continuum as a solid (dotted) thin 
line for T-ReCS (Spitzer) data. The /gj^com is estimated by the 
linear interpolation of this continuum flux at 9.7fj.m (in Fig. 3 
horizontal blue thick and red dot-dashed lines for T-ReCS and 
Spitzer data respectively). The error on the T9.7 was determined 
by the S/N of the continuum flux. More sophisticated methods 
as such spHne fit (Mason et al., 2012), PAHFIT (Smith et al., 
2007), or DecompIR (Mullaney et al., 2011) are also useful to 
measure T9 7 specially when the AGN is contaminated by extra- 
nuclear emission. However, as discussed in the previous section 
our sample are mostly AGN dominated. Moreover, these other 
methods require of wider wavelength ranges as they are based on 
a multiple component fitting and need a larger number of data- 
points to avoid degeneration of the results. However, no strong 
PAH features are detected in our sample (see Section 5). The op- 
tical depth has also been measured for the Spitzer spectra with 



SMART was developed by the IRS Team at Cornell University. 
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Fig. 4. Optical deptli T97 using SpitzerflRS and T-ReCS data. Filled 
symbols show the optical depth 7 using the linear interpolation in the 
8-13/im band and empty symbols show the optical depth T9.7 using the 
routine PAHFIT in the 5-35yum band (linked by short dotted lines). The 
size of the symbols scale with the strength of the PAHs. The continuous 
line shows the one-to-one relation. 



the same methodology. Note that PAHFIT may result more ap- 
propriate for Spitzer data given its longer wavelength range. We 
have however used the same method to measure T9.7 as for the 
T-ReCS spectra for consistency in the comparison of the silicate 

measurements. 

Strong PAHs could introduce uncertainties in the continuum 
determination using our linear interpolation (Spoon et al. , 2007). 
In order to check this, we have computed T9 7 using PAHFIT 
and the full wavelength range ( 5-35 iim) of the Spitzer/IRS 
spectra. Fig. 4 shows T9.7 for Spitzer spectra versus that of 
the T-ReCS spectra. Spitzer/IRS and T-ReCS data show sim- 
ilar T9 7 using the linear interpolation (see filled symbols in 
Fig. 4). The mean difference between T-ReCS and Spitzer T9.7 
is < T9 7(T - ReCS) - T9.7(Spitzer) > = 0.03. However, when 
we use a broader wavelength range and PAHFIT (see empty 
symbols in Fig. 4) we recover shallower optical depths for the 
Spitzer/IRS spectra (< T9.7(T - ReCS) - T9.7(Spitzer) > = 0.28). 
Thus, it seem that there is a systematic error of ~ 0.25 between 
the two methodologies. 

Only in two objects the value of T9.7 for the Spitzer/IRS 
spectra show a larger value than that found for T-ReCS, namely 
NGC 1808 and NGC 5135. Note that these two sources show the 
largest PAH features (represented with bigger symbols in Fig. 
4). In fact, we recover a lower value for both objects when using 
PAHFIT. In the following we will use the T9 7 obtained by the lin- 
ear interpolation for consistency with the T-ReCS data, although 
in Section 5.1 we discuss how T9.7 could affect the detection of 
the PAH feature. 

The apparent optical depths T9.7 measured are included in 
Table 3. Note that in NGC 3094 and NGC 7479 the bottom of 
the silicate absorption is very noisy so the determination of T9.7 



is poorly constrained. Thus, our measurement can be considered 
as lower limit. 

These T-ReCS data have already been presented in the lit- 
erature for a small number of individual objects. The optical 
depth at 9.7jum of NGC 3281 was already reported recently by 
Sales et al. (2011) using the same data. They found a value of 
T9.7 = 1.3 + 0.1, in goodagreementwithourresults. Diaz-Santos 
et al. (2010) also computed the silicate absorption of three ob- 
jects using T-ReCS data finding similar results. 

Most of the Spitzer/IRS spectra give similar T9 7 values to 
those measured from the T-ReCS data although the former tend 
to be smaller. When using PAHFIT and broader wavelength 
range these discrepancies become larger. This might reflect dif- 
ferent dust configurations (see Section 5.1). Large discrepan- 
cies are found only for NGC 1808, NGC 5 135, IC4518W and 
NGC 7479. For NGC 1808 and NGC 5135 the T9.7 estimates are 
lower than the Spitzer/IRS results. These two sources show the 
most contaminated Spitzer spectra of the sample (strong PAH 
features). NGC 1808 shows a lower value when using PAHFIT 
while NGC 5 135 shows an even larger value. 

We observed a wide range of silicate absorptions (t9.7=[- 
0.12,3.78]) with a mean value and a standard deviation of < 
T97 >= 1.20 ± 1.05. Only two sources show T9.7 < indicat- 
ing silicate features in emission. Among them the silicate fea- 
ture seen in emission for NGC 1365 is seen in absorption by 
Spitzer/IRS. The comparison with the literature is rather dif- 
ficult due to the different criteria to determine the continuum. 
However, we have compared T9.7 for Spitzer with the \2pm sam- 
ple reported by Wu et al. (2009), finding a good agreement for 
the 8 objects in common. 

The most significant features in the T-ReCS spectra are [SI V] 
at 10.5jum, the PAH features at 8.7 and 11.3yum, and [Nell] 
at 12.8jum''*. We have measured the equivalent width (EW) of 
[SIV] and PAH at 11.3jum dividing the integrated flux of the 
features by the adopted continuum near the feature ("feature" 
routine within IDL). We have computed the upper limit of the 
EWs accounting for the error in the continuum. We have also 
obtained EWs for these lines using Spitzer spectra. These EWs 
are also included in Table 3. The continuum luminosity at 12/im, 
vLy{l2pm), obtained for T-ReCS and Spitzer is also included in 
Table 3. 

5. Contributors to the MIR spectra of AGN 

5.1. The 11.3nm PAH feature 

The presence of PAHs has been associated to regions of star- 
formation, given the global correlations between star formation 
activity and PAH strength (e.g. Diamond-Stanic & Rieke, 2010; 
Roche et al., 1991; Smith et al., 2007; Wu et al., 2009). This 
emission is thought to originate in photo-dissociation regions 
(PDRs) where aromatic molecules are shielded from the radi- 
ation field produced by hot stars. 

It has been proposed that the hard radiation field produced 
by the AGN may destroy the PAH molecules (Genzel et al., 
1998; Siebenmorgen et al., 2004), as indicated by the decrease 
of the EW(PAH 1 1.3;um) with infrared colors (Wu et al., 2009). 
Nonetheless, Diamond-Stanic & Rieke (2010) found ll.3fim 
emission in most of their AGN in a sample of 35 Seyferts ob- 
served with Spitzer at scales of 1 kpc, finding suppression of the 

The short wavelength range prevent us to confirm the presence of 
the PAH feature at 8.7//m. Moreover, we have not attempted to measure 
the [Nell] line because in most cases it is too close to the limit of the 
wavelength range. 
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PAH features only for the 6.2, 7.7, and 8.6 fxm PAHs (see also 
Diamond-Stanic & Rieke, 2012). 

A strong PAH feature at 11 .S/jm is present only in NGC 1 808 
and NGC 7130 in the T-ReCS spectra. A hint of the PAH fea- 
ture is also present in other 4 AGN (NGC 1386, NGC 5643, 
NGC 7172, and NGC 7582). Thus, for 14 out of 20 sources their 
high spatial resolution T-ReCS spectra are free of signs of star- 
formation in the MIR, as traced by the PAHs. 

From visual inspection of the 2D T-ReCS spectra of these 
sources, the PAH feature is spatially extended for NGC 1808 and 
coincident with the nuclear emission for NGC 7130. Moreover, 
the central UV continuum of NGC 7130 consists of several knots 
of star-formation distributed in an asymmetric ring within 1 arc- 
sec, not spatially resolved by our data (Gonzalez Delgado et al., 
1998). Thus, the PAH features detected in NGC 7130 could be 
associated to this ring. NGC 1808 shows a circumnuclear star- 
burst region (e.g. Collison et al., 1994) together with an AGN 
(Veron-Cetty & Veron, 1985). The 3.3 yum PAH emission of 
NGC1808 was studied by Tacconi-Garman et al. (2005) using 
high spatial resolution images obtained with the ISAAC (VLT) 
instrument. They found PAH emission coincident with the nu- 
cleus supporting these results. 

Wu et al. (2009) described a method to estimate the AGN 
contribution when the PAH feature at 11.3jum is present. It is 
based on the constant ratio between the flux of the PAH feature 
and the continuum at 12;um iFi2^m^s f - '^'^■Ti^F pah) where the 
fraction of the AGN at MIR is Fagn = {Fn^im-Fntim,SF)IF\i^im. 
We have computed this value for the two objects with the 
strongest PAH emissions finding 87% and 93% for NGC 1808 
and NGC 7130, respectively. Therefore, the AGN contribution 
seems to be high even for the objects with the strongest PAH fea- 
tures. However, this method shows a 25% of error (see Section 
5.3 for further analysis on the AGN contribution of these two 
sources). 

Goulding et al. (2012), using SpitzerflRS spectra, found that 
star-formation dominates in most Compton-thick AGN. They 
argued that this might be an intrinsic characteristic of heavily 
absorbed AGN. Our sample includes 8 Compton-thick AGN 
but only NGC 7130 shows PAH emission in the T-ReCS spec- 
tra. However, as explained above, the T-ReCS spectrum of 
NGC 7 130 includes an unresolved ring of star-formation. Thus, 
it seems that the nuclear star-formation among Compton-thick 
AGN might be as uncommon as for the rest of the AGN. 
However, this needs to be confirmed for a larger sample of 
Compton-thick AGN. 

The lack of PAHs in the nuclear spectral spectra of AGN was 
already claimed by Siebenmorgen et al. (2004) using TIMMI2 
spectra. Honig et al. (20 10a) also reported a lack of PAH features 
in their high-spatially resolved VISIR spectra, detecting some 
minor emission of PAH at 11.3/iw only in 5 cases (in the two 
objects in common we also detect this tentative PAH emission). 

When we compare the T-ReCS and the Spitzer/IRS spectra, 
for those objects with strong PAH detected with T-ReCS, (i.e. 
NGC 1808 and NGC 7 130) the Spitzer PAH 11.3/zm feature is 
indeed present. Besides the PAH feature is detected in 16 other 
objects. The Spitzer spectra are free of PAH emission lines in 
only 4 objects, named NGC 3281, ESQ 103-G35, IC5063, and 
3C445. It seems that in most cases the PAH emission comes 
from emitting regions outside of 7- 1 30 pc (depending on the dis- 
tance) to the AGN. Considering that the Spitzer spectra have a 
spatial scale of 3 arcsec, that extension corresponds to anything 
between 0.6-0.8 arcsec in T-ReCS (depending on observation) 
and the 3 arcsec used for Spitzer spectra. This is consistent with 
previous studies that suggested that AGN lack strong PAH fea- 



tures (Genzel et al., 1998; Heman-Caballero et al., 2009; Honig 
et al., 2010a). 

However, note that the PAH at 11.3/im is embedded in the 
silicate absorption feature and thus in cases of high extinction 
the PAH feature could be attenuated (e.g. Brandl et al., 2006; 
Pereira-Santaella et al., 2010a). According to Smith et al. (2007), 
the optical depth at 1 1 .3/im is related to that at 9.1pm as ti 1 3^,,, = 
0.44 X Tcjjfim- For optical depths of rgj^^ = [1,2,3] (common 
optical depths in our sample, see Section 6) the intensity of the 
line could be reduced by a factor of / = x [0.64,0.41,0.26]. 
Thus, the attenuation of the PAH emission feature is expected 
to be large in some cases. Since the optical depths measured 
in our sample are similar for both T-ReCS and Spitzer spec- 
tra, we would expect that the PAHs must be attenuated in both. 
However, using PAHFIT for Spitzer spectra (see Section 4.2), 
the T9.7 optical depths are shallower in the Spitzer spectra than 
in T-ReCS. Thus, the lack of PAH features in the latter could be 
the result of the extra attenuation suffered by the T-ReCS spectra. 
This attenuation can be expressed as I - !„ x exp{-QAA{Tg,-i + 
AT9.7)), where At9,7 = T9,7(T - ReCS) - T9.7(Spitzer) is the dif- 
ference between the attenuation measured using T-ReCS and 
Spitzer spectra. Using the values of T9 7 given in Table 3, in most 
of the cases this extra attenuation cannot explain the lack of or 
weakness of the \\3pm PAH seen in T-ReCS data. In particu- 
lar, the upper limit on flie EW(PAH at II 3pm) for T-ReCS data 
could be explain as a result of the extra attenuation only in two 
galaxies, namely NGC 1386 and NGC 7479. Thus, the attenua- 
tion of the spectrum might not be responsible for the lack of PAH 
features in AGN. 

This lack of PAH emission in the proximity of the AGN 
(< 100 pc) indicates that the AGN emission is able to de- 
stroy the PAH molecules (Mason et al., 2007; Wu et al., 2009). 
Alternatively, the BH growth can continue after the gas fuel is 
not enough to maintain the star-formation activity (see the nu- 
merical simulations give by Hopkins, 2012). These delays have 
also been detected observationally (e.g. Davies et al., 2007). 
Therefore, the lack of PAH emission might be the result of a 
delay between the onset of the star formation rate and the black 
hole (BH) accretion rate. 

It is not easy to disentangle which of these two possibilities 
is the main reason for the lack of PAHs in AGN. If the PAHs 
are destroyed by the AGN we would expect that PAHs may sur- 
vive in the vicinity of low luminosity AGN. In support of this, 
we detected PAH features in the nuclear region of NGC 1808 
which is the object in our sample with the lowest X-ray lumi- 
nosity. Honig et al. (2010a) found that most AGN show lower 
continuum in the VlSlR than in the Spitzer data when the PAH 
are suppressed. They argued that this is in favor of the lack of 
star-formation as the main reason for the suppression of PAHs 
because the continuum of the MIR spectrum is roughly constant 
between VISIR and SpitzerflRS if the PAH destruction is an ef- 
fect of grain destruction. In our sample, although the continuum 
is roughly constant for some of the objects (e.g. Centaurus A), 
in some other objects the continuum suppression in the T-ReCS 
data is clear (e.g. NGC 7582). 

5.2. The [S IV] emission line at 10.5pm 

Another widely used AGN diagnostic line is the [SIV] at 10.5jum 
since this line is thought to be originated in the Narrow Line 
Region (NLR) (Dasyra et al., 2011, and references therein). 
However, the [SIV] at 10.5pm line arises from ions with ion- 
ization potential of 35 eV, similar to the optical line [OIII] at 
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Fig. 5. Absorption corrected 2-10 keV luminosity versus the llfiin lu- 
minosity, both in logarithmic units. Filled circles are T-ReCS data and 
open circles are Spitzer data. The dotted lines link Spitzer to T-ReCS 
data for the same object. Yellow stars and red circles mark the objects 
with the strongest [SIV] 10.5pm and PAH at ll.3pm emissions in the T- 
ReCS spectra, respectively. The blue dotted line shows the best fit found 
by Levenson et al. (2009) for AGN. Black dotted-dashed and solid lines 
represent the best fit found by Asmus et al. (201 1) for Starburst galaxies 
and AGN, respectively. 



5007A (see Trouille et al., 201 1). Thus, this line can be produced 
in star-forming regions as well as in AGN (Pereira-Santaella et 
al., 2010b). 

The [SIV] line is detected in 11 AGN and a hint of this 
line is also present in two other sources in the T-ReCS data 
(NGC3281, and NGC7582). In all these cases it is also de- 
tected with Spitzer. Interestingly, Spitzer detected this line in 
other five additional sources (NGC 1365, NGC 1808, Circinus, 
3C445, and NGC7479)'^ The lack of this line in these later 
objects could be due to the low S/N ratio of the T-ReCS spec- 
tra. In fact, they are the faintest emission lines, according to 
their EW([SIV]) in the Spitzer spectra. Besides, only NGC 1365 
shows a limit inconsistent with the EW reported with Spitzer 
data. Thus we claim that the origin of the [SIV] line is consistent 
with being nuclear in our sample (see also Pereira-Santaella et 
al., 2011). 

5.3. X-ray to MIR relation 

The mid-IR emission in AGN is originated in the reprocessing 
of the soft X-ray and UV radiation by dust located in the torus. 
At the same time, escaping high-energy photons could destroy 
PAH molecules at sub-kpc scales (Voit, 1992). This is supported 
by the lack of PAH emission in the T-ReCS spectra (see Section 
5.1). The reprocessed dust emission should then produce a sig- 
nificant correlation between X-ray luminosity (after correction 



'5 Note that this line is also detected in NGC 5728 but the T-ReCS 
spectra show a very low S/N (see Section 4.2). 



from X-ray absorption) and MIR luminosity. This correlation 
has already being found by other authors (Gandhi et al., 2009; 
Honig et al., 2010a; Horst et al., 2008; Ichikawa et al., 2012; 
Levenson et al., 2009). 

Fig. 5 shows the absorption corrected 2-10 keV luminos- 
ity L2-10 keV versus the ll/im continuum luminosity Lntim- X- 
ray errors represent the expected X-ray variability among AGN 
of 0.5 dex while the MIR errors are inferred from the photo- 
metric errors including 10% error for the total flux calibration. 
Spitzer luminosities (empty symbols) agree very well with the 
T-ReCS luminosities (filled symbols) for most of the sources 

(log(Li2,T-ReCs/Ll2,Spitzer) < 0.4). 

Most of the objects follow the same correlation found for 
AGN by Levenson et al. (2009) and Asmus et al. (2011) (blue- 
dotted and black-solid lines, respectively). Moreover, those ob- 
jects with significantly higher MIR emission compared to X-rays 
are also far from the correlation found by Asmus et al. (20 II) for 
starburst galaxies (dot-dashed line in Fig. 5). Note that the X-ray 
luminosity of NGC 7479 is not corrected for absorption although 
the object is Compton-thick. Thus, its X-ray luminosity could 
be up to 60 times higher than that reported here (see Panessa 
et al., 2006), in better agreement with the correlation for AGN. 
Interestingly, NGC 1808 (the closest to the starburst correlation) 
shows the strongest PAH emission in our T-ReCS sample, indi- 
cating that the AGN component of this source is not dominating 
the MIR emission. The locus of this object in the correlation in- 
dicates that the fraction of star-formation must be higher than the 
estimated (87%, see Section 5. 1 . In fact, the existence of an AGN 
in the nucleus of NGC 1808 has been under debate (e.g. Awaki 
et al., 1996; PhiUips, 1993). Jimenez-Bailon et al. (2005) stud- 
ied the X-ray emission of NGC 1808 with XMM-Newton and 
Chandra showing the co-existence of thermal diff'use plasma and 
non-nuclear unresolved point-like sources associated with the 
starburst activity, together with a low luminosity AGN. However, 
it is clear from our correlation that the AGN in NGC 1808 is not 
dominating the MIR emission. 

Thus, despite the heterogeneous sample and the large num- 
ber of highly obscured objects, the 12 jjin luminosity and the X- 
ray luminosities are still well correlated. Thus, except for the nu- 
cleus in NGC 1808, the AGN dominates the continuum emission 
at MIR frequencies. 



6. The silicate absorption/emission feature 

Wu et al. (2009) based in Spitzer/dsdSi found that the silicate fea- 
ture in Type-1 Seyferts is rather weak while in Type-2 Seyferts 
it is more likely to display strong silicate absorption. 

In the T-ReCS spectra, only the spectra of NGC 1365 and 
3C 445 show silicate features in emission while NGC 308 1 
shows very little absorption (see appendix A and Tgj in Table 
3). The other sources (i.e. 17 out of 20) show very strong sil- 
icate absorptions. Note that NGC 1365 and 3C445, the two 
objects showing the silicate feature in emission, are Type 1.5 
AGN. All the Compton-thick AGN in our sample ubiquitously 
show strong silicate absorption feature. This is in contrast with 
Goulding et al. (2012) who found that most of the Compton- 
thick AGN in their sample do not show such strong silicate ab- 
sorption. However, since our sample it is not complete, the statis- 
tics may well be biased by the selection of the individual ob- 
serving proposals of the targets in our sample. In particular none 
of the Compton-thick AGN with low optical depth reported by 
Goulding et al. (2012) have been observed by T-ReCS while we 
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Fig. 6. The optical depth at 9.1^m, Tg -j versus the logarithmic of the 
hydrogen column density at X-rays, Nh- Nh in units of cin^-. T-ReCS 
data are shown as blue-filled circles and the Spitzer/IKS data included 
in our work are shown as empty circles. The object without reported 
X-ray measurement is labeled as a short line in the Y-axis. Low limit on 
the Nh are shown as right arrows. Compton-thick sources reported by 
Goulding et al. (2012) are shown as open squares and AGN observed 
with VISIR/VLT by Honig et al. (2010a) are shown as black-filled stars. 
Constant gas-to-dust ratio (i.e. the uniform dusty torus) assuming the 
viewing angles of ; = 0, 40, 80° are shown with the continuous, dashed, 
and dot-dashed lines, respectively. 



6. 1 . Uniform dusty torus 

On the basis of a uniform (smooth), heavily-obscuring, sub par- 
sec scale torus (Granato & Danese, 1994; Granato et al., 1997; 
Pier & Krolik, 1992, 1993), the gas/dust-rich torus is responsi- 
ble for both the MIR and the X-ray gas absorption through the 
gas-to-dust relation. The constant gas-to-dust ratio can be esti- 
mated assuming Nh/Av - 1.93 x 10^'cm"^ and rgj - 0.07 x tx 
(Draine & Li, 2007), where tx is the total optical depth at X- 
rays. The solid line in Fig. 6 shows this constant value. For very 
small values of Nh it predicts no absorption, close to the actual 
values. Note however that some objects do show the silicate fea- 
ture in emission, which cannot be explained under these models. 
Nonetheless, as akeady shown by Goulding et al. (2012), this 
simple model is far from the measurements for most of the ob- 
jects. 

While the MIR emission is integrated over the whole dusty 
torus, the X-ray absorption is due only to the amount of ab- 
sorbing column density along the LOS. The influence of the Nh 
within the LOS can be estimated assuming a density profile at 
the viewing angle / (angle between the LOS and the plane of 
the torus). Nenkova et al. (2008a) assumed a Gaussian distribu- 
tion e"^, where cr is the width of the toroidal distribution. The 
constant gas-to-dust ratio for / - [0,40,80]° degrees is shown 
for cr = 45° (solid, dashed, and dotted-dashed black lines in 
Fig. 6, respectively). Only two of the sources could be explained 
by this relation (NGC7172 and NGC5506), assuming a highly 
inclined torus. Thus the observed Tgj feature is much shallower 
than predicted by this constant gas-to-dust relation. This was al- 
ready mentioned by Georgantopoulos et al. (201 1), showing that 
although ~70% of the sources with high optical absorption at 
9.7/im are classified as Compton-thick AGN, the opposite is not 
true, since many Compton-thick sources do not show significant 
absorption at 9.1 (see Fig. 6). 



do share three Compton-thick sources with strong silicate ab- 
sorption features'^. 

The X-ray hydrogen column density (Nh) of 20 of the 
sources reported in the literature are shown in Table 1 Fig. 6 
shows Tg 7 versus Nh in logarithmic scale. For comparison pur- 
poses we have included the high angular resolution data for AGN 
observed with VISIR/VLT (Honig et al., 2010a) as solid stars as 
well as the Compton-thick AGN (empty squares) observed with 
SpitzerflRS from the sample of Goulding et al. (2012). Note that 
filled symbols are chosen to emphasize high spatial resolution 
data (i.e. T-ReCS and VISIR). 

If we focus on the high spatial resolution data there is a trend 
in the sense of increasing rg 7 for higher Nh. This relation has al- 
ready been seen by several authors (e.g. Alonso-Herrero et al., 
2011; Honig et al., 2010a; Shi et al., 2006; Wu et al., 2009). 
Moreover, for log(NH) > 23 the silicate absorption shows a wide 
range of values from low or no absorption to deep silicate ab- 
sorptions in agreement with the results presented by Wu et al. 
(2009). In the next sections we will try to model this relation- 
ship under two assumptions: (1) the uniform dusty torus model 
and (2) the clumpy torus models. 



Note that we have five sources in common with Goulding et al. 
(2012) although two of them did not have high S/N T-ReCS data for the 
spectral analysis. 
" Note that there ai'e not X-ray data published for NGC 3094. 



6.2. Clumpy torus 

The clumpy torus models of Nenkova et al. (2008a,b) predict the 
silicate feature in emission even for high values of the optical 
depths of the individual clouds (xy > 100). This effect can also 
explain the low silicate absorption observed for some Compton- 
thick sources in our sample. In order to test this possibility we 
have used the libraries of CLUMPY"* models, a code modeling 
the AGN clumpy dust torus emission (Nenkova et al., 2008a,b). 
We have obtained a set of spectral energy distributions (SEDs) 
for cr = 45°, Y - RoiitlRinner - 200 (where Ri„„er and are the 
inner and outer radii of the torus, respectively), an exponential 
slope of the radial distribution of clouds q - 2 (see Nenkova et 
al., 2008b, for more details), ry = 5 - 150, and the number of 
clouds along the equator of the torus N,, - 2 - 25 clouds. We 
chose Y to use the lowest maximum temperature and therefore 
the maximum silicate depth. Hereinafter we assume cr - 45°. 
Alonso-Herrero et al. (2011) has measured cr - [30 - 60]° for 
a sample of Type-2 Seyferts using a Bayesian method to fit the 
CLUMPY models (see also Asensio Ramos & Ramos Almeida, 
2009). Note that this parameter is not crucial since the silicate 
absorption is insensitive to the angular size of the torus when the 
LOS incidence angle of the ray is smaller or close to the size of 
the torus (i.e. Type-2 Seyferts, see Nenkova et al., 2008b). Under 
this approximation, the total optical depth produced by the torus 
is the sum of the optical depths of the individual clouds along 
the LOS (tx - N XTy, where is the number of clouds along 



http://www.pa.uky.edu/clumpy 
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the LOS) and the number of clouds A^, similarly to the classical 
torus, is 



N = N„e' 



(2) 



We have applied the same technique described in Section 4.2 
to obtain the silicate feature strength for these SEDs. Fig. 7 (dot- 
dashed black lines) shows T9 7 for these models for an inclination 
; = (i.e. along the plane of the torus) and No of 5 and 25 clouds. 

Nonetheless, this number of clouds is still an estimate 
of the mean number of clouds along the LOS. Statistically for 
a viewing angle / and a mean equatorial number of clouds No 
the ray might intercept a range of clouds that might differ sub- 
stantially from the average. We have estimated the range of 
following the prescription given by Nenkova et al. (2008a) (see 
their Appendix). Using Poissonian distributions the probability 
that a ray intercepts a number of clouds k normalized to the prob- 
ability to intercept the mean value of clouds is: 



Pn 



m 
J] 



(3) 



the range of the number of clouds (k) is obtained as the minimum 
and maximum for which the probability is smaller than 10%. 
The dotted lines toward the higher (lower) Nh in Fig. 7 repre- 
sents the expectation for the largest (lowest) number of clouds 
with a probability above 10%. 

Many of the Compton-thick AGN are well within the bound- 
aries of this model if we are observing them close to the plane 
of the torus and the number of clouds is high (A^o ~ 5 - 25). The 
colored thick lines show how this model changes with the incli- 
nation angles for the two No- The model reproduces well all the 
mildly obscured X-ray objects (with Nh < lO^^cm"^) assuming 
a mean number of clouds of A^o > 5. This is in good agreement 
with the parameters of the Bayesian fit to CLUMPY models re- 
ported by Ramos Almeida et al. (2009), Ramos Almeida et al. 
(201 1), and Alonso-Herrero et al. (201 1). 

However, the model fails to reproduce a dozen sources with 
high value for T97. Nenkova et al. (2008a) showed that irre- 
spective of the optical depth, the absorption feature produced 
by a clumpy torus is never deep (rg 7 < 1). This is confirmed 
by other clumpy models developed (e.g. Honig & Kishimoto, 
2010b). Alonso-Herrero et al. (2011) also noted that it is possi- 
ble that for the Seyfert galaxies with the deepest silicate features 
a clumpy medium in a torus-like configuration may not be appro- 
priate to explain the observations. Levenson et al. (2007) showed 
that deep silicate features require that the source is embedded in 
dust that is both optically and geometrically thick. They showed 
that a geometrically thick shell structure is able to produce T9.7 
as high as 4. One possibility then, is that these objects show dust 
screen geometries in addition to the clouds of the clumpy torus. 

6.3. Contribution from extended dust to tlie observed silicate 
feature 

On the basis of the uniform dust torus models produced from ra- 
diative transfer theory (Dullemond & van Bemmel, 2005; Fritz 
et al., 2006; Schartmann et al., 2005), the MIR spectrum of ob- 
scured AGN included in our sample are expected to be domi- 
nated by a power-law and significant silicate absorption features 
at A ~ 9.7yum. However, as discussed in the previous sections, 
many AGN in our sample show silicate absorption features shal- 
lower than predicted for the homogeneous dusty model and of- 
ten deeper than predicted for the clumpy model. Several authors 
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Fig. 7. The optical depth at 9.7/ira, Tgj versus the logarithmic of the 
hydrogen column density at X-rays, Nh- Nh in units of cm"-. Symbols 
as in Fig. 6. Dot-dashed lines are the CLUMPY model predictions for 
I = 0° and No = [5, 25] clouds. Dotted lines show how the curves evolve 
for the minimum and maximum number of clouds with a probability 
higher than 10% of intercepting the LOS. The colored thick line show 
how the curve evolve when the viewing angle changes from ; = 0° 
(black) until ; = 80° (yellow). The dotted lines are the model including 
the Poisson probability to detect a number of clouds along the ray (see 
text). 



have claimed that the dust seen at MIR frequencies might not 
be co-spatial with the torus gas measured at X-rays (Goulding et 
al., 2012; Sturm et al., 2005). In particular, the inclination angle 
of the host galaxy, dust lanes, or galaxy minor or major merg- 
ers may have significant eff'ects on the AGN emission at MIR 
(Alonso-HeiTero et al., 2011; Deo et al., 2007, 2009; Goulding 
& Alexander, 2009; Goulding et al., 2012; Malkan et al., 1998; 
Matt et al., 2000). Although not many indications of extended 
contributors to the emission can be seen in our spectra, most of 
the sources in our sample are extended (see Section 4.1). In this 
section we discuss the possible contribution of the host galaxy 
to the optical depth at 9.1 ^m. 

To study a possible extra-nuclear origin for the silicate ab- 
sorption feature we have looked for the superb spatial resolution 
that HST can achieve. Broad band HST images are available in 
the MAST archive''' for 12 AGN in our sample (10 with the 
F606W filter and 2 with the F814W filter). Narrow-band filters 
have been used in 5 cases whenever broad-band filter images 
were not found. For the remaining 5 objects without HST im- 
ages, we have used the available catalog of images in NED^" 
(NGC3094, NGC4418, Circinus, IC4518W, and 3C445). We 
have compiled the axial ratio of the galaxies as the b/a ratio, 
where b and a are the minor and mayor axis respectively, to 
measure the inclination of the host galaxy (see Table 1). Axial 
ratios of b/a < 0.5 are considered as edge-on galaxies. Dust 
lanes or extended dust across the nucleus are present in 8 ob- 



" http://archive.stsci.edu/ 
http://ned.ipac.caltech.edu 
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Fig. 8. The optical depth at 9Jfjm, Tgj versus the logarithmic of the hy- 
drogen column density at X-rays, Nh- Nh in units of cm^^. Spitzer/IKS 
data of our sample are removed here for clarity. Edge-on galaxies are 
shown as big red squares, galaxies with dust-lanes are shown with big 
yellow stars, and mergers are shown as big green circles. Note that edge- 
on galaxies are shown only when dust-lanes are not present for clarity 
of the plot. 



jects and the AGN is hosted in edge-on galaxies in 9 objects. In 
10 objects no particular inclination or dust features are found. 
We have also checked for the morphological properties of the 
host galaxies of the sample. Five galaxies are in mergers, namely 
NGC4945 (Aalto et al., 1991), Circinus (Aalto et al., 1991), 
IC4518W (Vorontsov-Velyaminov, 1977), NGC7130 (Aalto et 
al., 1991), and NGC7479 (Laine & Heller, 1999; Martin et al., 
2000). These properties are included in Table 1. Note that for 
comparison purposes the same properties hava been looked for 
in the Goulding et al. (2012) and Honig et al. (2010a) samples. 

Fig. 8 shows Tgj versus Nh where we mark those AGN 
hosted in edge-on galaxies as red squares, AGN with nuclear 
dust-lanes as yellow stars, and mergers as green circles. For clar- 
ity we show only edge-on galaxies without nuclear dust-lanes. 
We note that all the objects with high Tgj are hosted in galax- 
ies meeting one or more of these three properties. In general, 
clumpy torus models can only explain face-on galaxies free of 
either dust lanes and/or participating in mergers, although this 
analysis cannot rule out the host galaxy as the responsible for 
the observed Tgj, even for low Tg j. Only three objects show one 
of these properties without a strong (rg 7 < 0.4) silicate absorp- 
tion features (IC4329A, NGC3227, and NGC424). Therefore, 
we conclude that even at the high angular resolutions ( 100 pc) of 
our observations, dust in the host galaxy can play a role on the 
overall dust extinction seen at MIR frequencies, and particularly 
it does when deep silicate absorption features are observed. 

7. Conclusions 

Here we present an analysis of 20 high spatial resolution MIR 
AGN spectra observed with T-ReCS (Gemini observatory) at 



scales of ~100 pc. The T-ReCS data have been processed with 
the pipeline RedCan. RedCan has been developed by our group 
for the data reduction of imaging and spectroscopic data from 
the CanariCam instrument, recently commissioned at the GTC 
(La Palma, Spain). This data set has been complemented with 
Spitzer spectra probing kpc scales as well as X-ray information 
for the sources. The main results of the study are: 

1. On the AGN contribution to the MIR spectrum at scales of 
-100 pc. We find that the AGN dominates the MIR emis- 
sion as seen in the high-spatial resolution T-ReCS spectra 
presented here. This is in agreement with the relatively faint 
ll.3fim PAH emission observed in most of the AGN studied 
here. 

2. On the origin of the silicate feature. The observed strength 
of the nuclear (a few hundred pc) 9.7/im optical depth could 
be well fitted with clumpy models with a No ~ 5 - 25 num- 
ber of clouds but only for sources with little of mild silicate 
absorption. However, sources with the deepest 9.1iim optical 
depths cannot be reproduced with clumpy models. In these 
cases we have found that the host galaxies contribute to the 
nuclear component since they are detected in objects show- 
ing optical dust lanes and/or they are involved in a merger 
process. 

3. On the importance of the spatial resolution. Intense \ \.3^m 
PAH feature emission is detected in most of the Spitzer[[RS> 
spectra. This indicates that, at least in our sample of AGN, 
most of the star formation activity, as probed by this tracer is 
taking place in regions at least 100 pc away from the AGN. 
Interestingly, in most of the AGN the 12jUot Spitzer contin- 
uum luminosity match those observed by T-ReCS. 
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Appendix A: Catalogue of spectra 







Fig. A.l. T-ReCS spectra of the sample are shown as open circles. The SpitzerflRS spectra are the red lines. The continuum fits to the T-ReCS and 
IRS spectra are shown as solid and dotted lines, respectively. Positions of the most prominent emission lines are shown by dashed lines. 
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Fig. A.l. T-ReCS spectra of te sample. Continued. 



